Strong tunneling and charge quantization in S-I-N Coulomb blockade structures 
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We study the charge of a smaU normal island connected by a tunnel junction to a superconducting 
lead. Unlike the N-I-N case, the steps of the Coulomb staircase remain sharp even if the conductance 
of the tunneling barrier exceeds e^ /h. One can observe the transition from sharp steps to the smeared 
ones by applying magnetic field to destroy the superconductivity. 
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The charge of a small conducting island weakly cou- 
pled to an electron reservoir (lead) by a tunneling bar- 
rier shows a step-like dependence on the gate voltage due 
to the well-known phenomenon of Coulomb blockade Jd It 
was shown in Refs. that the quantization of charge 
is not precise: the processes of virtual tunneling between 
the lead and the grain result in smearing of charge steps. 
The smearing is especially dramatic in the regime of 
strong tunneling, when the conductance G of the tun- 
neling barrier is large, G ^ e^/h. The evolution of the 
charge quantization with the increasing G was studied 
in semiconductor quantum dots, where the junction ci 
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ductance can be tuned between G = and G e /hi 
Charge quantization was also observecfl in a metallic 
"electron box" device, yielding a sharp Coulomb stair- 
case at small G. In a metallic device, one can achieve 
the tunneling juncticp, conductance G e^//i, Ref. ||. 
In this case, theorjitltil predicts only weak oscillations of 
charge, with amplitude ~ e exp(— /iG/2e^). However, the 
layout of a metallic device does not allow one to vary the 
junction conductance, and the crossover between weak 
and strong tunneling cannot be observed in one device. 

In this paper we study quantum fluctuations of charge 
of a normal grain connected by tunnel junction to a su- 
perconducting lead. We show that due to the existence 
of the gap A in the lead the steps of the Coulomb stair- 
case remain vertical at any G. The fluctuations result 
in a finite slope of the plateaus of the staircase. Un- 
like the N-I-N case, the slope may remain small even 
at G » e'^/h. One can suppress the gap by apply- 
ing a magnetic field|3 thus achieving a crossover between 
the limits of weak and strong charge fluctuations in one 
sample. The main result of the paper is illustrated by 
Fig. 0, obtained under the assumption A ^ Ec, where 
Ec = e^/2G is the charging energy of the electron box, 
and C is its total capacitance. 

To account for the Coulomb interactions, we include 
the charging energy term (Q — q)^ /2G into the Hamil- 
tonian H of the system. The external charge q here is 
proportional to the gate voltage. Then by statistical av- 
eraging of the obvious relation dH/dq = {q — Q)/G, we 
find the average charge of the grain 



Here F{q) = — T In Z{q) is the free energy. 

To find the partition function Z{q), we use the stan- 
dard effectiiie action approach by Ambegaokar, Eckern 
and Schon£3 In this formalism an auxiliary Hubbard- 
Stratonovich field <p(t) is introduced to replace the 
quadratic in Q interaction term by a linear one. The 
electronic degrees of freedom are then traced out, result- 
ing in the following expression for the partition function: 
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with 



Zm ^ Dip{t) exp < - / -rrrdT - St[ 
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Here the summation over winding numbers m accounts 
for the discreteness of chargeji3 and /? is inverse temper- 
ature. The tunneling contribution St[ip\ to the effective 




Q{q) = {Q)=q-C 



dF{q) 
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FIG. 1. Coulomb staircase for a normal grain connected 
to a superconducting lead. The steps remain vertical at any 
conductance. The plot corresponds to the case A 3> Ec, and 
G = {e^/h)A/Ec, see Eq. (0). Note that G > e'^/h. 
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action was evaluated for a wide junction in Rcf. |2[ It is 
non-local in time r; however, at large A ^ Ec its contri- 
bution reduces to a mere renormalization of capacitance 
C in the action in Eq. (^). The magnitude of the renor- 
malization was found in Ref. In the case of N-I-S 
junction the renormalized capacitance is 



(7 = C + ?iG/2A. 
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In this approximation the action in Eq. is quadratic, 
and we can easily find the dependence of Zm on the wind- 
ing number: 



Zm = Zq CXp 
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To study Z{q) at low temperatures T ^ e^/C*, it is con- 
venient to transform the partition function (||) with the 
help of the Poisson summation formula. Equations (|^) 
and d) then yield 



Z{q) - Zo 



27rCT 



+ 00 

E 



exp 



{ne — cif 



2CT 



(6) 



The largest term in sum (^) corresponds to an integer n 
which is closest to q/e. At T this term dominates 
the sum, and using Eqs. (||) and (^, we find 



Q(g)=elnt( ^ + i 
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where Int (a;) is the integer part of x. This is the central 
result of the paper; see also Fig. ^ At G — + the second 
term in Eq. (Q) vanishes, leading to a perfect staircase. 
At non-zero G the second term describes the finite slope 
of the plateaus. According to Eq. (0) the steps in Q{q) 
remain sharp at any conductance, in contrast to the case 
of a normal lead. Even at G ^ jfi the heights of the 
steps e/(l + EcfiG / /S.e^) are still close to their nominal 
value e. 

The fact that Q{q) is discontinuous at half-integer q/e 
is due to the density of states in the superconductor van- 
ishing at the Fermi level, and therefore does nat rely on 
the assumption A ^ Ec- Indeed, it is knowrO that the 
problem of charge fluctuations in an electron box can be 
mapped onto an effective spin-i Kondo problem. It fol- 
lows from Eq. (30) of Ref. ^ that the height of the charge 
step is 5Q = 2e/i, where ^ is the renormalized value of 
the spin of the Kondo impurity. In the case of a normal 
metal the spin is completely screened, = 0, and the 
steps are smeared completely, 5Q — Q. If at least one of 
the leads is a superconductor, the spin /j, is not comletely 
screened, and the step height 5Q remains finite. 

At finite but low temperature T <^ e^/C, the steps 
are smeared slightly. To account for the smearing, one 
needs to retain two terms, with n — Int(g/e) and n = 
Int^q/e) + 1, in the sum (|^). In this approximation the 
charge becomes 
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One can easily check that in the limit of zero temperature 
this result reproduces Eq. (Q). 

The effective action technique used here is applica- 
ble for tunnel junctions of wide area, where the con- 
ductance is distributed over a large number A'^ ^ 1 of 
transverse modes. At finite N one may have to take into 
account the possibility of two-electron (Andreev) tun- 
neling. Although the steps in the dependence Q{q) re- 
main vertical, the plateau conductance due to such pro- 
cesses acquires a correction, which can be estimated as 
SQ ^ {hG/e'^yN~^. The diffusion of electrons inside 
the grain may enhance this correction. The result cau 
be expressed in terms of effective number of channelsLEI 
Nqs ^ N, which still remains very large due to the small- 
ness of the Fermi wavelength in metals. 
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